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D | An updated analysis using about 1.5 million events recorded at y/s = Mz with 

the DELPHI detector in 1994 is presented. Eighteen infrared and collinear safe 
event shape observables are measured as a function of the polar angle of the 
thrust axis. The data are compared to theoretical calculations in 0(a 2 ) includ- 
ing the event orientation. A combined fit of a s and of the renormalization scale 
in 0(a 2 ) yields an excellent description of the high statistics data. 

The weighted average from 18 observables including quark mass effects and cor- 
relations is a s (M|) = 0.1174±0.0026. The final result, derived from the jet cone 
energy fraction, the observable with the smallest theoretical and experimental 
uncertainty, is 

ot B (M|) = 0. 1180 ± 0.0006(exp.) ± 0.0013(hadr.) ± 0.0008(scale) ± 0.0007(mass). 

Further studies include an a s determination using theoretical predictions in the 
next-to-leading log approximation (NLLA), matched NLLA and 0(a 2 ) pre- 
dictions as well as theoretically motivated optimized scale setting methods. 
The influence of higher order contributions was also investigated by using the 
method of Pade approximants. Average a s values derived from the different 
approaches are in good agreement. 
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1 Introduction 

This paper presents a highly improved test of second order perturbation theory and 
an improved measurement of a s (M|). It is based on progress in next-to-leading order 
QCD calculations of oriented event shape distributions |T|]. Furthermore, the DELPHI 
data used in this analysis are much improved in both their statistical and systematic 
precision compared with those of previous DELPHI publications |2|||. The distributions 
of 18 different infrared and collinear safe hadronic event observables are determined from 
1.4 million hadronic Z-decays at various values of the polar angle $t of the thrust axis. 
The *&t dependence of all detector properties are taken into account, thus achieving the 
best possible experimental precision. 

The precise experimental data are fully consistent with the expectation from second 
order QCD. A two parameter fit to each of the distributions measured at different polar 
angles $t allows an experimental optimization of the 0(a 2 ) renormalization scale giving 
a consistent set of eighteen a s (M^) values. For most of the distributions the largest 
uncertainty on the a s (M|) values is due to hadronization corrections and not to renor- 
malization scale errors. Any artificial increase Q of the uncertainty of a s (M|) due to a 
large variation of the renormalization scale is avoided so that the degree of precision to 
which QCD can be tested remains transparent. An average value of a s (M|) is derived 
taking account of the correlations between the values obtained from the 18 distributions. 

A number of additional studies have been performed to check the reliability of the 
a s (M|) results obtained from experimentally optimized scales. In one of these studies 
'optimized' renormalization scales as discussed in the literature are used to determine 
a s (M|) in second order pertubation theory. The different methods applied for choosing 
an optimized scale lead to consistent results for the average value of a s (M|). However, 
the scatter among a s (M|) values from the individual distributions is smaller for the 
experimentally optimized scales than that obtained using theoretically motivated scale 
evaluation methods. The correlation between the renormalization scales obtained with 
the different methods is also investigated. 

Further determinations of a s (M|) are performed by using all orders resummed cal- 
culations in the next-to- leading logarithmic approximation (NLLA). Here two different 
methods are applied. In the first case the pure NLLA predictions are confronted with 
the data in a limited fit range. In the second method a s is determined using matched 
NLLA and 0(a 2 ) calculations. For both methods the renormalization scale is chosen to 
be [i = M z . Both methods lead to average a s values consistent with the average value 
obtained in 0(a 2 ) with experimentally optimized renormalization scales. The agreement 
between the results of the pure NLLA fits and those of the 0{a 2 ) is emphasized. A closer 
inspection of the fits in matched NLLA and 0(a 2 ) to the very precise data reveals a so 
far unreported problem with this method in that the trend of the data deviates system- 
atically from the expectation of the matched theory. 

The selection of hadronic events and the correction procedures applied to the data are 
described in Section ||. Section § introduces the investigated event shapes and compares 
the expectations from various fragmentation models. Section ^ contains the comparison 
with angular dependent second order QCD and a detailed discussion of the determina- 
tion of a s (M%). Section [5] summarises determinations of a s (M|) using the renormalization 
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scales discussed in the literature. Section |6| discusses results obtained by applying Pade 
approximants for the extrapolation of the pertubative predictions to higher orders. Sec- 
tion |7| discusses results from applying NLLA. The heavy quark mass correction of the 
a s (M<l) values derived from experimentally optimized renormalization scales is described 
in Section ^. The final results are summarized in the last section. 

2 Detector and Data Analysis 

In this analysis the final data measured with the DELPHI detector in 1994 at a centre- 
of-mass energy of a/s = Mz are used. The statistics of the 1994 data is fully sufficient for 
the accurate QCD studies described in this paper. DELPHI is a hermetic detector with a 
solenoidal magnetic field of 1.2 T. The detector and its performance have been described 
in detail in ||. The following components are relevant to this analysis: 

• the Vertex Detector, VD, measuring charged particle track coordinates in the plane 
perpendicular to the beam with three layers of silicon micro-strip detectors at radii 
between 6.3 and 11 cm and covering polar angles, with respect to the e~ beam 
between 37° and 143°; 

• the Inner Detector, ID, a cylindrical jet chamber with a polar angle coverage from 
17° to 163°; 

• the Time Projection Chamber, TPC, the principal tracking detector of DELPHI, 
which has 6 sector plates in both the forward and backward hemispheres, each with 
16 pad rows and 192 sense wires, inner and outer radii of 30 cm and 122 cm and 
covers polar angles from 20° to 160°; 

• the Outer Detector, OD, a five layer drift chamber at 198 to 206 cm radius covering 
polar angles between 43° and 137°; 

• two sets of forward planar drift chambers, FCA and FCB, with 6 and 12 layers 
respectively and overall polar angle coverages of 11° to 35° and 145° to 169°; 

• the High Density Projection Chamber, HPC, a lead-glass electromagnetic calorime- 
ter with a very good spatial resolution located inside the DELPHI coil between 208 
cm and 260 cm radius and covering polar angles between 43° and 137°; 

• the Forward Electromagnetic Calorimeter, FEMC, comprising two lead-glass arrays, 
one in each endcap, each consisting of 4500 lead-glass blocks with a projective ge- 
ometry, and covering polar angles from 10° to 36.5° and from 143.5° to 170°; 

• The hadron calorimeter, HAC, an iron-gas hadronic calorimeter outside the coil, 
consisting of 19 to 20 layers of streamer tubes and 5 cm thick iron plates also used 
as flux return, whose overall angular coverage is from 11.2° to 168.8°. 

2.1 Event Selection 

Only charged particles in hadronic events were used. They were required to pass the 
following selection criteria: 

• momentum, p, greater than 0.4 GeV/c, 

• Ap/p less than 100%, 

• measured track length greater than 30 cm, 

• track polar angle between 16° and 164°, 

• impact parameter with respect to the nominal interaction point within 4 cm per- 
pendicular to and 10 cm along the beam. 
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Hadronic events were selected by requiring: 

• at least 5 charged particles, 

• the total energy of charged particles greater than 12%-y/i where the pion mass has 
been assumed for all particles, 

• the charged energy in each hemisphere of the detector, defined by the plane perpen- 
dicular to the beam, Ehemis, greater than 3%y/s, 

• the polar angle of the thrust axis between 90.0° and 16.3°. 

In total about 1.4 million events satisfy these cuts. The selection efficiency is 92%. 
Since the thrust axis does not distinguish between forward and backward directions, it 
is chosen such that cos^t > 0. is called the event orientation. The data are binned 
according to the event orientation into eight equal bins of cos$y between 0.0 and 0.96. 
With the exception of the eighth bin, the thrust axis is well contained within the detector 
acceptance. 

2.2 Correction Procedure 

The contamination of beam gas events, 77 events and leptonic events other than t + t~, 
is expected to be less than 0.1 % and has been neglected. The influence of r + r~ events 
which have a pronounced 2-jet topology and contain high momentum particles has been 
determined by the KORALZ model || treated by the full simulation of the DELPHI 
detector DELSIM |7j and the standard data reconstruction chain. The t + t~ contribu- 
tions have been subtracted from the measured data according to the relative rate of t + t~ 
(0.46% ± 0.03%) and hadronic events. 

The observed data distributions were corrected for kinematic cuts, limited acceptance 
and resolution of the detector as well as effects due to reinteractions of particles inside 
the detector material. The simulated data were processed in the same way as the real 
data. The correction for initial state photon radiation has been determined using events 
generated by JETSET 7.3 PS || with and without initial state radiation as predicted by 
DYMU3 ||. For any given observable Y the bin- by-bin correction factor C(Y, cos^t) is 
calculated as: 



Particles with a lifetime larger than 1 ns were considered as stable in the generated 
distributions. The bin widths were chosen on the basis of the estimated experimental 
resolution so as to minimize bin-to-bin migration effects. 

For the evaluation of systematic errors the cuts for the track and event selections 
were varied over a wide range, including additional cuts on the momentum imbalance, 
etc. Variation of the tracking efficiency has been considered by discarding 2% of the 

1 Thrust and the thrust axis have been calculated applying the algorithm included in the JETSET package 13. 

2 For the comparison of event shape observables with QCD predictions in all orders resummed next-to-leading-log 
approximation, the distributions have been integrated over Differing from the event selection criteria listed above, the 
hadronic events were selected if the polar angle of the thrust axis satisfied 40.0° < i?t < 90.0° for these angular integrated 
distributions. 
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accepted tracks at random. The influence of uncertainties in the momentum resolution 
was estimated by applying an additional Gaussian smearing of the inverse momenta 
of the simulated tracks. From the stability of the measured distributions a systematic 
uncertainty has been computed as the variance with respect to the central value. As 
the systematic error is expected to grow in proportion to the deviation of the overall 
correction factor from unity, an additional relative systematic uncertainty of 10% of this 
deviation has been added quadratically to the above value. 

3 Measured Event Shapes and Comparison with 
Fragmentation Models 

This analysis includes all commonly used shape observables which have a perturbative 
expansion known at least to next to leading order. 

3.1 Definition of the Observables 



Thrust T is defined by [10 



T = max^^ , (2) 

where pi is the momentum vector of particle i, and is the thrust axis to be determined. 

Major M and Minor m are defined similarly, replacing Ht in the expression above by 
the Major axis nMaj, which maximizes the momentum sum transverse to ut or the Minor 
axis nMin = ^Maj x Ut respectively. 



The oblateness O is then defined by ITT 



= M-m . (3) 
The C-parameter C is derived from the eigenvalues A of the infrared-safe linear mo- 



mentum tensor 4J 12 



Qid = I V ^ (A) 

£*Ip*I V &l 

C = 3 • (AiA 2 + A 2 A 3 + A 3 Ai) . (5) 
Here p\ denotes the z-component of Pk- 

Events can be divided into two hemispheres, a and b, by a plane perpendicular to the 
thrust axis n^. With M a and Mb denoting the invariant masses of the two hemispheres, 
the normalized heavy jet mass pn, light jet mass p^, the sum of the jet masses ps and 
their difference po can be defined as 



max(M a 2 ,M b 2 y 

PH = ^2 ( 6 ) 

vis 

min(M Q 2 ,M b 2 ) 

PL = 7^2 {') 



vis 
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PS = PH+ PL 



PD = Ph - PL 



where 



(9) 



(10) 



and the energy of the particles % has been calculated assuming pion mass for charged and 
zero mass for neutral particles. 



Jet broadening measures have been proposed in |TJ| . In both hemispheres a and b the 
transverse momenta of the particles are summed thus: 

_ E, ga , fe \PiXn T \ 

B a ,b ~ „ ^ f-ZTj • (11) 

2 Ei \Pi\ 

The wide jet broadening B max , the narrow jet broadening B min , and the total jet 
broadening B sum are then defined by 



B m ax = max(5 a , B b 



(12) 



B„ 



mm(B a , B b ) 



(13) 



B, 



B m ax B n 



The first order prediction in perturbative QCD vanishes for both /?l and B 
fore these observables cannot be used for the determination of a s . 



(14) 
mm- There- 



Jet rates are commonly obtained using iterative clustering algorithms |14| in which a 
distance criterion or a metric j/y, such as the scaled invariant mass, is computed for all 
pairs of particles % and j. The pair with the smallest y^ is combined into a pseudoparticle 
(cluster) according to one of several recombination schemes. The clustering procedure is 
repeated until all of the yij are greater than a given threshold, the jet resolution parameter 
y cut . The jet multiplicity of the event is defined as the number of clusters remaining; the 
n-jet rate R n (y C ut) is the fraction of events classified as n-jet, and the differential two-jet 
rate is defined as 



D 2 (y, 



R2(y cu t) - R2(y cu t - Ay 



cut I 



cut ) 



Ay, 



(15) 



cut 



Several algorithms have been proposed differing from each other in their definition of 
y^ and their recombination procedure. We apply the E0, P, P0, JADE |15" |, Durham 
T6| , Geneva |14[ and the Cambridge algorithms ||17| . The definitions of the metrics y^ 



and the recombination schemes for the different algorithms are given below. 



In the E0 algorithm y^ is defined as the square of the scaled invariant mass of the pair 
of particles i and j: 
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VV = —ff2 • ( 16 ) 



(Pi ±Pj) 2 

VIS 

The recombination is denned by: 

E k 

E k = Ei + Ej , p k = (pi + pj) , (17) 

\Pi + Pj | 

where Ei and Ej are the energies and pi and pj are the momenta of the particles. 
In the P algorithm yij is defined by Eq. ([16]), and the recombination is defined by 



Pk = Pi+Pj , E k = \p k \ . (18) 

The PO algorithm is defined similarly to the P algorithm, however the total energy 
E vis (Eq. [[(J) is recalculated at each iteration for the remaining pseudoparticles. 

In the JADE algorithm, the definition of is 

2E l E J (l- cos e tj ) 

vis 

where 9ij is the angle between the pair of particles i and j. 

For the Durham algorithm is given by 

2min(£? £ 2 )(i_ cos fl ) 
Vij = (20) 



2 

vis 



and for the Geneva algorithm by 



8^(l-cos%) 

For the algorithms given by Equations ([19]), (p0|) and ( pT]) the recombination is done 
by adding the particles four-momenta. 

The recently proposed Cambridge algorithm [0] introduces an ordering of the particles 
i and j according to their opening angle, using the ordering variable 

^• = 2(l-cos%) (22) 

and is defined by Eq. (|20|) . The algorithm starts clustering from a table of N^j 
primary objects, which are the particles' four-momenta, and proceeds as follows: 

1. If only one object remains, store this as a jet and stop. 

2. Select the pair of objects i and j that have the minimal value of the ordering variable 
Vij and calculate y^ for that pair. 

3. If < y cut then remove the objects i and j from the table and add the combined 
object with four-momentum Pi+Pj. If y^ > y cu t then store the object i or j with the 
smaller energy as a separated jet and remove it from the table. The higher energy 
object remains in the table. 

4. go to 1. 
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The energy-energy correlation EEC |18] is defined in terms of the angle Xij between 



two particles i and j in an hadronic event: 

^(x) = ^EE§^ / w-x*w > (23) 

a. AT ■ ■ J - J vis J 



N i,j 



x- 



where N is the total number of events, Ax is the angular bin width and the angle % is 
taken from \ = 0° to x = 180°. 

The asymmetry of the energy-energy correlation AEEC is defined as 

AEEC(x) = EEC(180°-x)-EEC( X ) • (24) 

The jet cone energy fraction JCEF [Tj| integrates the energy within a conical shell of 
an opening angle x about the thrust axis. It is defined as 

JCEF ^ = ^EEi- J . (25) 



Ax 



where Xi is the opening angle between a particle and the thrust axis vector Ht, whose 
direction is defined here to point from the heavy jet mass hemisphere to the light jet mass 
hemisphere. Although the JCEF is a particularly simple and excellent observable for the 
determination of a s , it has been rarely used until now in experimental measurements. 
Within an 0{a 2 s ) analysis, the region 90° < x < 180°, corresponding to the heavy jet 
mass hemisphere, can be used for the measurement of a s . The distribution of the JCEF 
is shown in Figure |l[ Hadronization corrections and detector corrections as well as the 
next-to-leading order perturbative corrections are small. This allows a specially wide fit 
range to be used. 



3.2 Fragmentation Models 

QCD based hadronization models, which describe well the distributions of the event 
shape observables in the hadronic final state of e + e~ annihilation, are commonly used for 
modelling the transition from the primary quarks to the hadronic final state. Perturba- 
tive QCD can describe only a part of this transition, the radiation of hard gluons and 
the evolution of a parton shower. For a determination of the strong coupling constant a s 
one has to take account of the so-called fragmentation or hadronization process, which 
is characterized by a small momentum transfer and hence a breakdown of perturbation 
theory. Several Monte Carlo models are in use to estimate the size of the hadroniza- 
tion effects and the corresponding uncertainty. The most frequently used fragmentation 
models, namely JETSET 7.3 PS {§, ARIADNE 4.06 and HERWIG 5.8c jn| have 



been extensively studied and tuned to DELPHI data and to identified particle spectra 
from all LEP experiments in As discussed in detail in [^] all models describe the 
data well. Examples of the measured hadron distributions are presented in Figures [I] 
and @. The increased systematic accuracy of the data is partially due to the fact that 
the t?T dependence of the detector corrections is explicitly taken into account. The $t 
dependence of the detector corrections is shown, for instance, in Figure H for two of the 
observables studied. Detailed tables of the individual event shape distributions including 
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their statistical and systematic errors will be made available in the HEPDATA database 



For figures of the distributions see also 



Before theoretical expressions describing parton distributions can be compared with 
experimental data, corrections have to be made for hadronization effects, i.e. effects 
resulting from the transition of the parton state into the observed hadronic state. For 
the global event shape observables this transition is performed by a matrix P, where Py 
is the probability that an event contributing to the bin j of the partonic distribution 
will contribute to the bin i in the hadronic distribution and is computed from a Monte 
Carlo model. This probability matrix has been applied to the distributions from 0(a 2 s ) 
perturbative theory D pert .(Y, cos^t) to obtain the distributions for the predictions of the 
observed final state Dh a dr.(Y, cosi?t): 

D ha dr.{y, COS$ T )i = ^ P ij( Y i COS $ T ) D perL (Y, COS0 T )i ■ (26) 

3 

In the case of the JCEF, EEC and AEEC, which are defined in terms of single particles 
and pairs of particles, respectively, bin-by-bin correction factors Cnadr. similar to that 
described above for the detector effects have been computed such as: 

D hadr XY, cos $ T )i = C Hadr ,(Y,cos'd T ) i D pert XY,cos'd T )i . (27) 

Our reference model for evaluating hadronization effects is the JETSET 7.3 Parton 
Shower (PS) Generator, which has been modified with respect to the heavy particle de- 
cays to obtain a better description of the heavy particle branching fraction. This modified 
version is denoted by JETSET 7.3 PS D in the following. The tuned parameters have 
been taken from |2"2]], where the updated tuning procedure is described in detail. 



In order to estimate the systematic error of the hadronization correction, the analy- 
sis was repeated using alternative Monte Carlo generators with different hadronization 
models. In addition, the parameters for the JETSET PS were varied. A description of 
the models and their differences can be found for example in |22| . The alternative models 
used are ARIADNE 4.06, HERWIG 5.8c as well as version 7.4 JETSET PS IB. All these 



models have been tuned to DELPHI data [22|. For our standard Monte Carlo program 



we applied also an alternative tuning to the DELPHI data which includes Bose-Einstein 
correlations, not included in the reference tuning. Whereas the number of hard gluons 
predicted by second order QCD matrix elements is simulated by the hadronization mod- 
els |[25|| , additional soft gluons are produced within the parton shower cascade, controlled 



by the JETSET PS parameter Q , which describes the parton virtuality at which the 
parton shower is stopped. To account for the sensitivity of the shape observables with 
respect to the additional soft gluons, Q has been varied from 0.5 GeV to 4.0 GeV. 



The systematic error of a s originating from hadronization corrections is then estimated 
as the variance of the fitted a s values obtained by using all the hadronization corrections 
mentioned above. Further studies have been made to investigate the influence of the 
main fragmentation parameters of the JETSET PS model by varying them within their 
experimental uncertainty. It has been found that this contribution to the uncertainty of 
a s in general is less than one per mille, and has been neglected. 




Figure 1: left part: Measured 1-T distribution integrated over cos^t- The upper part shows the detector correction including effects 
due to initial state radiation. The part below shows the size of the hadronization correction. The width of the band indicates the 
uncertainty of the correction. In the central part the measured 1-T distribution is compared to the expectation from four hadronization 
generators, JETSET 7.3 PS D with DELPHI modification of heavy particle decays, JETSET 7.4 PS, ARIADNE 4.06 and HERWIG 
5.8c. Also shown is the 1-T range used in the QCD fit. The lower part shows the ratio (Monte Carlo simulation-data) /data for the four 
hadronization generators. The width of the band indicates the size of the experimental errors, right part: Same curves as shown in the 
left part but for JCEF integrated over cos$t- ^ 




Figure 2: left part: Measured 1-T distribution in two bins of cos^t- The upper part shows the detector corrections in the two cos^t 
bins. The part below shows the size of the relative hadronization correction in the two cos^t bins with respect to the average correction. 
In the central part the measured 1-T distributions are compared to JETSET 7.3 PS D. The lower part shows the ratio (Monte Carlo 
simulation-data)/data for the two cos^t bins, right part: Same curves as shown in the left part but for JCEF in two bins of cos$t- 
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4 Comparison with Angular Dependent Second Or- 
der QCD using Experimentally Optimized Scales 

The evaluation of the 0(a 2 ) coefficients is performed by using EVENT2 a pro- 
gram for the integration of the 0(a 2 ) matrix elements. The algorithm is described in p|. 

Using this program, one can calculate the double differential cross-section for any 
infrared and collinear safe observable Y in e + e~ annihilation as a function of the event 
orientation: 



1 d 2 a(Y, cos-^r) , 2 \ a 

= a a (jr) ■ A(Y, cos0 T ) 



a tot dYd cos -&T 

i -2/ 2 

+ <w 



£(Y,costf T ) + f 2vr/3 ln(x M ) - 2) A(Y, cos $ T ) 



,(2* 



where a s = a s /27r and /3q = (33 — 2n/)/127r, n/ is the number of active quark flavours 
and Otot is the one loop corrected cross-section for the process e + e~ — > hadrons. The 
event orientation enters via i?t, which denotes the polar angle of the thrust axis with 
respect to the e + e~ beam direction. The renormalization scale factor x M is defined by 
/i 2 = x^Q 2 where Q = Mz is the centre-of-mass energy. A and B denote the 0(a s ) and 
0(a 2 ) QCD coefficients, respectively. Alternatively, the double differential cross-section 
can be normalized to the partial cross-section in each cos$t interval: 

R{Y,COS$t) = ~, ^ ~AVA Q ( 29 ) 

\ a cos vt J dY a cos vt 
which is more appropriate for the study of residual QCD effects. 

The strong coupling a s at the renormalization scale \x is in second order perturbative 
QCD expressed as 



fa lnln^ 



= 1 - -k^f , (30) 



(5) 

where A = A^ is the QCD scale parameter computed in the Modified Minimal Subtrac- 
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tion (MS) scheme for rif — 5 flavours and (3\ = (153 — 19n/)/247r 

The renormalization scale \i is a formally unphysical parameter and should not enter 
at all into an exact infinite order calculation ||27|| . Within the context of a truncated finite 
order perturbative expansion for any particular process under consideration, the defini- 
tion of /i depends on the renormalization scheme employed, and its value is in principle 
completely arbitrary. This renormalization scale problem has been discussed extensively 
in the literature p|p7|j2~8i . 

The traditional experimental approach to account for this problem has been to mea- 
sure all observables at the same fixed scale value, the so-called physical scale x M = 1 or 
equivalently ji 2 = Q 2 . The scale dependence has been taken into account by varying \i 
over some wide ad hoc range, quoting the resulting change in the QCD predictions as 
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theoretical uncertainty. 

However, the approach of choosing x^ = 1 has a severe disadvantage. If we consider 
the ratio of the 0(a s ) and the 0(a 2 s ) contributions to the cross-section, defined as 



B(Y) + A(Y)(2TTf3 Q Mx ti )-2) 
r NL o(Y) = — ] - , (31) 



A(Y) 

we find for many observables quite large values for the second order contributions. In 
some cases this ratio can have a magnitude approaching unity, indicating a poor conver- 
gence behavior of the 0(a 2 ) predictions in the MS scheme which would quite naturally 
result in a wide spread of the measured a s values. This has indeed been observed in 
previous analyses using 0(a 2 ) QCD [p|, |29| , |30 |. 



Several proposals have been made which resolve the problem by choosing optimized 
scales according to different theoretical prescriptions [0T],P^,f53|]. These methods have been 
discussed with some controversy P7| and until now no consensus has been achieved. The 
only approach for determining an optimized scale value, which does not rely on specific 
theoretical assumptions, is the experimental evaluation of an optimized 0(a 2 ) scale value 
x^ for each measured observable separately. This strategy has therefore been chosen to 
be the primary method. For previous analyses including experimentally optimized renor- 
malization scales see for example |29| , |3^| . The theoretical approaches for choosing an 
optimized renormalization scale value are studied in detail in Section |5[ As will be shown 
later, the approach of applying experimentally optimized scales yields an impressive con- 
sistency of the a s (M|) measurements from different observables. 

The procedure applied here, is a combined fit of a s and the scale parameter x M . In the 
past this strategy suffered from a poor sensitivity of the fit with respect to x M for most 
of the observables. Due to the high statistics and high precision data now available, one 
may expect a better sensitivity at least for some of the observables under consideration. 

We determined a s (M|) and the renormalization scale factor x M simultaneously by 
comparing the corrected distributions for each observable Y with the perturbative QCD 
calculations corrected for hadronization effects as described in the previous section. The 
theoretical predictions have been fitted to the measured distributions R(Y, cos fix) by 
minimizing % 2 , defined by using the sum of the squares of the statistical and systematic 
experimental errors, with respect to the variation of and x M . 

The fit range for the central analysis, i.e. including the experimental optimization of 
x M , was chosen according to the following considerations: 

• Requiring a detector acceptance larger than 80%, the last bin in cos^t was excluded 
in general, i.e. the fit range was restricted to the interval < cos^t < 0.84 which 
corresponds to the polar angle interval 32.9° < "$t < 90.0°. 

• Acceptance corrections were required to be below about 25% and the hadronization 
corrections to be below ~ 40%. 

• The contribution of the absolute value of the second order term r NLO (Y) as defined 
in Eq. ( |3TD was required to be less than onef] over the whole fit range. 



^This requirement restricts the fit interval only for the total jet broadening observable B sum , which yields rather large 
O(atg) contributions for any choice of the renormalization scale value. 
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• The requirement that the data can be well described by the theoretical prediction, 
i.e. x 2 / n df is approximately 1 and stable over the fit range. 

• Stability of the a s - measurement with respect to the variation of the fit range. 

For the analysis with a fixed renormalization scale value x^ — l, the requirement that 
X 2 / n df is approximately 1 can in general not be applied, since it would cause an unrea- 
sonably large reduction of the fit range for many observables. The thrust distribution 
for example could be fitted only over a range of at most three bins. For further details 



see also Section fO|. The fit ranges for this analysis as well as for the analyses applying 
theoretically motivated scale setting methods have therefore been chosen identical to the 
analysis with experimentally optimized scale values, regardless of the \ 2 values of the 
fits. 



4.1 Systematic and Statistical Uncertainties 

For each observable the uncertainties from the fit of a s (M§) and of have been de- 
termined by changing the parameters corresponding to a unit increase of x 2 . In the case 
of asymmetric errors the higher value was taken. 

The systematic experimental uncertainty was estimated by repeating the analysis with 
different selections to calculate the acceptance corrections as described in Section ||. Addi- 
tionally, an analysis was performed including neutral clusters measured with the hadronic 
and/or electromagnetic calorimeters. The overall uncertainty was taken as the variance 
of the individual a s (M|) measurements. 

An additional source of experimental uncertainty arises from the determination of the 
fit range, which has been estimated by varying the lower and the upper edge of the fit 
range by ±1 bin, respectively, while the other edge is kept fixed. Half of the maximum 
deviation in a s (M§) has been taken as the error due to the variation of the fit range and 
has been added in quadrature. 

The hadronization uncertainty was determined as described in Section [3[ 

The total uncertainty on a s (M§) is determined from the sum of the squares of the 
errors listed above. 



Uncertainties due to Missing Higher Order Calculations 

An additional source of theoretical uncertainty arises due to the missing higher order 
calculations of perturbative QCD. It is commonly assumed, that the size of these un- 
certainties can be estimated by varying the renormalization scale value applied for the 
determination of a s (M|) within some 'reasonable' range [j35 ]. The choice of a 'reasonable' 



range involves subjective jugdement and so far no common agreement about the size of 
this range has been achieved. Furthermore, this commonly used approach has been crit- 
icized in the literature M. According to [|J any artificial increase of the uncertainty of 
a s (M|) due to a large variation of the renormalization scale should be avoided so that 
the degree of precision to which QCD can be tested remains transparent. It should be 
pointed out that no such additional uncertainty is required to understand the scatter of 
the measurements from a large number of observables if experimentally optimized renor- 
malization scale values are applied. This will be demonstrated in the following section. 
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Other procedures for estimating uncertainties due to missing higher order corrections 
have been suggested, in particular the comparison of a s (M|) values obtained by applying 
different reasonable renormalization schemes or by replacing the missing higher order 
terms by their Pade Approximants [[JIJ. Both strategies have been studied. By compar- 
ing the size of the uncertainties derived applying these methods (see e.g. Table 11), a 
variation of x M between 0.5 • x e * p and 2 • x e * p seems justified to obtain an estimate of these 
uncertainties. Similar or identical ranges have for example also been chosen in PJ51] . 



4.2 Results 

The results of the fits to the 18 event shape distributions are summarized in tables [| 
and ^| and shown in Figure 0. It should be noted that for all observables the normalized 
X 2 is about one for a typically large number of degrees of freedom (n^/ = 16 — 236, see 
Table []]). The individual errors contributing to the total error on the value of a s are 
listed in Table |2|. Among the observables considered the JCEF yields the most precise 
result. 



For comparison, the data have also been fitted in 0(a 2 s ) applying a fixed renormaliza- 
tion scale value x^ = 1. The results of these fits are summarized in table ||] and shown in 
Figure |j. As can be seen from the x 2 l n df values of the fits, the choice of x M = 1 yields 
only a poor description of the data for most of the observables, for many observables the 



Observable 


Fit Range 


cos 


d T 


Range 






X 2 /n d f 


n d f 


EEC 


28.8° - 


151.2° 





- 





84 


0.0112 ±0 


0006 


1 


02 


236 


AEEC 


25.2° - 


64.8° 





- 





84 


0.0066 ± 


0018 





98 


75 


JCEF 


104.4° - 


169.2° 





- 





84 


0.0820 ± 


0046 


1 


05 


124 


1 - T 


0.05 - 


0.30 





- 





84 


0.0033 ± 


0002 


1 


24 


89 





0.24 - 


0.44 





- 





84 


2.30 ±0 


40 





90 


33 


c 


0.24 - 


0.72 





- 





84 


0.0068 ± 


0006 


1 


02 


82 


BMax 


0.10 - 


0.24 





- 





84 


0.0204 ± 


0090 





89 


47 


Bsum 


0.12 - 


0.24 





- 





84 


0.0092 ± 


0022 


1 


19 


40 


Ph 


0.03 - 


0.14 





- 





84 


0.0036 ± 


0004 





63 


54 


PS 


0.10 - 


0.30 





- 





36 


0.0027 ±0 


0019 





82 


16 


Pd 


0.05 - 


0.30 





- 





84 


2.21 ±0 


38 


1 


02 


68 


u 2 


0.07- 


0.25 





- 





84 


0.048 ± 


020 





85 


68 


d po 

u 2 


0.05 - 


0.18 





- 





84 


0.112 ±0 


048 


1 


02 


68 


u 2 


0.10 - 


0.25 





- 





84 


0.0044 ± 


0004 


1 


00 


47 


T» Jade 
u 2 


0.06 - 


0.25 





- 





84 


0.126 ±0 


049 


1 


05 


75 


T\ Durham 
u 2 


0.015 - 


0.16 





- 





84 


0.0126 ±0 


0015 





92 


96 


T\ Geneva 
u 2 


0.015 - 


0.03 





- 





84 


7.10 ±0 


28 





84 


19 


r\ Cambridge 
U 2 


0.011 - 


0.18 





- 





84 


0.066 ±0 


019 





98 


145 



Table 1: Observables used in the 0(a 2 ) QCD fits. For each of the observables the fit 
range, the range in cos^r, the measured renormalization scale factor x M together with the 
uncertainty as determined from the fit, the x 2 / n df and the number of degrees of freedom 
ndf are shown. In the case of asymmetric errors the higher value is given. 
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description is even unacceptable. 



More details concerning the QCD fits are presented in Figures |5| to || Figure ^| shows 
the values of a s (M|) and the corresponding values of A% 2 , i.e. the change of y 2 with 
respect to the optimal value, for the fits as a function of the scale lg(x^) for some of 
the investigated observables. The shape of the A% 2 curves indicates that for most dis- 
tributions the renormalization scale has to be fixed to a rather narrow range of values 
in order to be consistent with the data. For most of the observables the renormalization 
scale dependence of a s (M|) is significantly smaller in the region of the scale value for the 
minimum in x 2 / n df than for the region around x M = 1. It should however be noted, that 
even for observables exhibiting a strong scale dependence of a s (M|) , e.g. D^ eneva ^ the 
a s (M|) value for the experimentally optimized scale value is perfectly consistent with the 
average value. 

Figures || to |9| contain a direct comparison of the data measured at various bins in 
cos i?t with the results of the QCD fits. The measured dependence on both cos $t and 
the studied observable are precisely reproduced by the fits. 

At this point a comparison with the results from applying a fixed renormalization scale 



value Xu = 1 seems appropriate. Figure [10] shows QCD fits to the data with experimen- 



Observable a s (M|) Aa s (Exp.) Aa s (Hadr.) Aa s (Scale.) Aa s (Tot.) 



EEC 





1142 


± 





0007 


± 





0023 


± 





0014 


± 





0028 


AEEC 





1150 


± 





0037 


± 





0029 


± 





0100 


± 





0111 


JCEF 





1169 


± 





0006 


± 





0013 


± 





0008 


± 





0017 


1 -T 





1132 


± 





0009 


± 





0026 


± 





0023 


± 





0036 








1171 


± 





0028 


± 





0030 


± 





0038 


± 





0056 


c 





1153 


± 





0021 


± 





0023 


± 





0017 


± 





0036 


Btvlax 





1215 


± 





0022 


± 





0031 


± 





0013 


± 





0041 


Bsum 





1138 


± 





0030 


± 





0032 


± 





0030 


± 





0053 


Ph 





1215 


± 





0014 


± 





0029 


± 





0050 


± 





0060 


PS 





1161 


± 





0014 


± 





0018 


± 





0016 


± 





0033 


Pd 





1172 


± 





0013 


± 





0034 


± 





0007 


± 





0038 


D E0 

u 2 





1165 


± 





0027 


± 





0029 


± 





0017 


± 





0044 







1210 


± 





0018 


± 





0026 


± 





0009 


± 





0033 


Dl 





1187 


± 





0019 


± 





0021 


± 





0036 


± 





0046 


T) Jade 
u 2 





1169 


± 





0011 


± 





0020 


± 





0028 


± 





0040 


T\ Durham 
u 2 





1169 


± 





0013 


± 





0016 


± 





0015 


± 





0026 


T\ Geneva 
u 2 





1178 


± 





0052 


± 





0075 


± 





0295 


± 





0309 


-■-.Cambridge 
U 2 





1164 


± 





0008 


± 





0023 


± 





0004 


± 





0025 



Table 2: Individual sources of errors of the a s (M|) measurement. For each observ- 
able, the value of a s (M|) , the experimental uncertainty (statistical and systematic), the 
uncertainty resulting from hadronization corrections, the theoretical uncertainty due to 
scale variation around the central value x e ® p in the range 0.5 • x e * p < x^ < 2 • x e ® p and 
the total uncertainty are shown. 
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tally optimized and with fixed renormalization scale values x M = 1 for the observables 
1—T,ph and ■ It can be seen that the slope of the experimental distributions cannot 
be described by the theoretical prediction applying x M = 1. The same observation has 
also been made for other observables. Good description of the data can in general only 
be achieved within the small kinematical region where the fit curve intersects with the 
data. Demanding x 2 / n df — 1 the thrust distribution for example can only be described 
within a maximum range of three bins in 1-T. It should be noted, that in contrast to 
fits applying experimentally optimized scales, the stability of a s (M|) with respect to the 
choice of the fit range is in general quite poor. This observation is new and due to the 
fact that the data used in this analysis have both smaller statistical and systematical 
errors than in previous publications. (See also Section [Of ) and |39]. 



Combining the 18 individual results from the a s (M|) measurements applying experi- 
mentally optimized renormalization scale values by using an unweighted average yields 



a s (M|) = 0.1170 ± 0.0025 

whereas the corresponding average for the measurements using the fixed scales x^ = 1 
is a s (Mf) = 0.1234 ± 0.0154. For the experimentally optimized scales the scatter of the 
individual measurements is significantly reduced. 

The consistency of the individual measurements using experimentally optimized scales 
is clearly shown by the good x 2 / n df — 9.6/17 for the unweighted average. This value is 



Observable ajM£) Aa s (Scale) Aa s (Tot.) xyndf 



EEC 





1297 


± 





0037 


± 





0042 


10.7 


AEEC 





1088 


± 





0015 


± 





0050 


2.04 


JCEF 





1191 


± 





0012 


± 





0024 


7.7 


1-T 





1334 


± 





0042 


± 





0051 


25.9 


O 





1211 


± 





0065 


± 





0077 


2.38 


c 





1352 


± 





0043 


± 





0053 


12.0 


BMax 





1311 


± 





0073 


± 





0083 


1.67 


Bsum 





1403 


± 





0056 


± 





0071 


8.1 


Ph 





1325 


± 





0036 


± 





0049 


5.1 


PS 





1441 


± 





0055 


± 





0062 


2.16 


Pd 





1181 


± 





0012 


± 





0039 


1.54 


D E0 
u 2 





1267 


± 





0033 


± 





0052 


1.35 


d po 

u 2 





1265 


± 





0026 


± 





0041 


1.31 


D p 

u 2 





1154 


± 





0019 


± 





0036 


5.35 


~p\ Jade 
u 2 





1249 


± 





0030 


± 





0042 


1.53 


j^Durham 





1222 


± 





0034 


± 





0046 


3.47 


r~\ Geneva 
u 2 





0735 


± 





0071 


± 





0116 


120. 


-pvCambridge 
U 2 





1202 


± 





0021 


± 





0033 


1.32 



Table 3: Results of the a s (M|) measurements using a fixed renormalization scale — 1. 
For each observable, the value of a s (M^) , the uncertainty from the variation of the scale 
between 0.5 < x^ < 2, the total uncertainty and the x 2 / n df of the fit are shown. 
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computed on the basis of the total uncertainty (experimental and hadronization uncer- 
tainty) without considering an additional renormalization scale error. For the average of 
the fixed scale x^ — 1 measurements x 2 / n df = 168/17, thus in this case the individual 
measurements are clearly inconsistent with each other. This inconsistency may be under- 
stood to arise from the imperfect description of the data for the fits with x M = 1, which 
cause a s , the parameter of the fit, not to be well defined. 

The idea behind the common analysis of such a large number of observables is to 
optimize the use of the information contained in the complex structure of multi-hadron 
events. Errors due to the corrections for hadronization effects may be expected to can- 
cel to some extent in the averaging procedure. To test this expectation the analysis of 
each of the individual 18 observables is repeated by performing hadronization corrections 
with all hadronization generators described in Section 3.2. This results in 7 times 18 
individual a s values. As a first test for each of the 18 observables the unweighted av- 
erage value of a s from the seven models is evaluated. The average value of the 18 a s 
values is a a (M|) = 0.1177 ± 0.0029. In a second step for each of the 7 hadronization 
models an unweighted average of the corresponding 18 a s values is calculated. Finally an 
unweighted average of the 7 average values for the different hadronization models is com- 
puted resulting in a s (Ml) = 0.1177 ± 0.0016. The result confirms that the scatter of the 
average values due to different assumptions for hadronization corrections is significantly 
smaller than the uncertainty of ±0.0025 of the mean value from 18 individual observables. 

The correlations between the a s values obtained from the different observables must be 
taken into account in order to calculate their weighted average. Since the correlations are 
mostly unknown, the exact correlation pattern cannot be worked out reliably. Therefore, 
we use a recently proposed method which makes use of a robust estimation of the 
covariance matrix and has been used for example in [37] . Here it is assumed that different 



measurements i and j are correlated with a fixed fraction p e g of the maximum possible 
correlation C™ ax : 

= pesC™ ax i^j, with c;p' a.a, . (32) 

For p cS = the measurements are treated as uncorrelated, for p eff = 1 as 100% 
correlated entities. When x 2 < n df the measurements are assumed to be correlated and 
the value p e s then is adjusted such that the x 2 is equal to the number of degrees of 
freedom 

X 2 (pcs) = - x)(xj - x)(C~ l )ij = n-l = n df . (33) 

When x 2 > n df it is assumed that the errors of the measurement are underestimated 
and will therefore be scaled until x 2 — n df is satisfied. 

Applying this method to the 18 observables studied, the weighted average (see also 
Fig. |3|) yields: 

a s (Mf) = 0.1168 ± 0.0026 

with p e fj being 0.635. Both the central value and its uncertainty are almost identical to 
the unweighted average and the r.m.s. quoted above, which in itself is a remarkable result. 
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It should be noted, that the method applied for the calculation of the weighted average 
does not necessarily lead to the smallest possible error. In [57 it is shown for example, 



that the error of the weighted average is increased if less significant measurements are 
included. Within this analysis however, the a s measurements from all individual ob- 
servables have been considered, regardless of their significance. This is motivated by 
the fact that the errors of the a s measurements quoted in Table |2] and in the following 
subsections contain all uncertainties which can be evaluated from a careful experimental 
analysis. However, the spread of the a s measurements may not be explainable by the 
individual uncertainties alone, there may be additional uncertainties, which cannot be 
derived from a single observable. Therefore the above averaging procedure has been ap- 
plied and robustness of the error estimate has been preferred instead of minimizing the 
error. Still this error estimate may not cover a possible general shift of the measured 
average with respect to the true a s (M|) value. Apart from a better theoretical under- 
standing, today such a shift could only be inferred by comparing to different types of 
calculation, like resummed or Pade approximation, which are presented in Sections |6| and 

4.3 Additional Cross Checks 

An underlying assumption for the 0(a 2 s ) QCD fits to the shape observables is that 
the value of a s is approximately independent of a variation of the renormalization scale 
within the fit range. To check this assumption, a cross check for the differential two jet 



rate observables has been performed following a suggestion in ||38|| . For these observables 
the QCD fits have been repeated, allowing the renormalization scale to vary proportion- 
ally to y m ti i-e. /i 2 = x^Dcut^fs. The differences in the a s (M|) determination have been 
found to be of the order of a few per mille for the individual jet rate observables and to 
be less than two per mille for the average of these observables. 

A further investigation has been performed for all observables of Table |l|. The fit 
range listed in Table Q has been divided into two separate, approximately symmetrical 
regions, allowing a maximum overlap of one bin. a s (M|) has been determined applying 
experimentally optimized scales for both regions independently. The fits were successful 
for all observables except £)^ eneva y where the resulting fit ranges were too small to allow 
the fits to converge. Good agreement of the two a s values measured for each observ- 
able is found. In a further step, the two a s (M|) values have been combined for each 
observable according to their statistical weight. These a s values have been combined 
by calculating a weighted average as described before. The resulting average value of 
a s (M§) = 0.1168 ± 0.0025 is identical to the value determined from the standard pro- 
cedure. No systematical trend of the two values of the renormalization scales found for 
the two fit ranges (dominated by two respectively three jet events) is observed. Further 
information on this study can be found in |39] . 
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a s (M z 2) = 0.1168 ± 0.0026 

X 2 /n df =6.2/17 

Peff = 0.635 



0.06 0.08 0.1 0.12 0.14 0.16 0.18 

a s (M z 2) 



Figure 3: Results of the QCD fits applying experimentally optimized scales for 18 event 
shape distributions. The error bars indicated by the solid lines are the quadratic sum 
of the experimental and the hadronization uncertainty. The error bars indicated by 
the dotted lines include also the additional uncertainty due to the variation of the 
renormalization scale due to scale variation around the central value x e * p in the range 
0.5 • x e ® p < x^ < 2 ■ x e * p . Also shown is the correlated weighted average (see text). The 
X 2 -value is given before readjusting according to Eq. [J3L 
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Figure 4: Results of the QCD fits applying a fixed renormalization scale x^ — 1 . The 
error bars indicated by the solid lines are the quadratic sum of the experimental and 
the hadronization uncertainty. The error bars indicated by the dotted lines include also 
the additional uncertainty due to the variation of the renormalization scale around the 
central value x e * p from 0.5 • x e ® p < x^ < 2 ■ x e ® p . Also shown is the correlated weighted 
average. It has been calculated assuming the same effective correlation p eff = 0.635 as 
for the fit results applying experimentally optimized scales. The x 2 / n df for the weighted 
average is 71/17, where the \ 2 given corresponds to the value before adjusting p cS . In 
order to yield x 2 / n df = 1, the errors have to be scaled by a factor f err = 3.38. 




Figure 6: (left side) QCD fits to the measured thrust distribution for two bins in cos-#t- (right side) Measured thrust distribution at 
various fixed values of 1 — T as a function of cos$t- The solid lines represent the QCD fit. 




Figure 7: Same as Figure 6 but for the jet cone energy fraction JCEF 
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Figure 8: Same as Figure 6 but for the energy energy correlation EEC. 




Figure 9: Same as Figure 6 but for the differential two jet rate with the Jade Algorithm D^ 6 . 



to 

CI 




Figure 10: Comparison between O(a^) QCD fits with x^ — 1 and experimentally optimized renormalization scales for the observables 
1 — T, p H and D% '. DELPHI data and theoretical predictions are shown averaged over cos-# T - The full (dotted) line correspond to fits 
with experimentally optimized renormalization scale values x e * v [x^ = 1), respectively. The theoretical prediction is averaged over the 
bin width and drawn through the bin center. 
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5 Scale Setting Methods from Theory 

Several methods for the choice of an optimized value for the renormalization scale 
have been suggested by theory. See e.g. [|28| for an overview. In this section we compare 
the a s (M|) measurements, using renormalization scales predicted by three different ap- 
proaches with the a s (M|) measurements using experimentally optimized scales: 

(i) The principle of minimal sensitivity (PMS) : Since all order predictions should be 



independent of the renormalization scale, Stevenson |31| suggests choosing the scale to 



be least sensitive with respect to its variation, i.e. from the solution of 

da 



dx^ 



. (34) 



(ii) The method of effective charges (ECH) : The basic idea of this approach |32| 
is to choose the renormalization scheme in such a way that the relation between the 
physical quantity and the coupling is the simplest possible one. In 0(a 2 s ) , where the 
ECH approach is equivalent to the method of fastest apparent convergence (FAC) |32| 
the scale is chosen in such a way that the second order term in Eq. ( p8|) vanishes: 

B(Y, cos # T ) + (2n/3 la(x fl ) - 2)A(Y,cos$ T ) = . (35) 

(Hi) The method of Brodsky, Lepage and MacKenzie (BLM) ||33|| : This method follows 
basic ideas in QED, where the renormalized electric charge is fully given by the vacuum 
polarization due to charged fermion-antifermion pairs [p8[| . In QCD it is suggested to 
fix the scale with the requirement that all the effects of quark pairs be absorbed in the 
definition of the renormalized coupling itself. In 0(a 2 s ) this amounts to the requirement 
that x^ is chosen in such a way that the flavour dependence nf of the second order term 
in Eq. |28| is removed: 







dn 



f 



. (36) 



B(Yn f ) + (2n(3 \n(x fl )-2)A(Y) 

) n f =5 

The results of the a s (M|) measurements for the individual observables applying the 
different scale setting prescriptions are listed in table [|. The weighted averages for the 
different methods yield: 

(i) PMS method : 

a s (M 2 ) = 0.1154 ± 0.0045 (x 2 /n df = 19/17) 

(ii) ECH method : 

a s (M 2 z ) = 0.1155 ± 0.0044 (x 2 /n df = 19/17) 

(in) BLM method : 

a s (M|) = 0.1174 ± 0.0068 (x 2 /n d f = 29/13) 

to be compared with a s (M|) = 0.1168 ± 0.0026 ( x 2 / n df = 6.2/17) using the experimen- 
tally optimized scales. 
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The weighted averages for the different theoretical methods are in agreement with the 
average using experimentally optimized scales. The scatter of the individual measure- 
ments is lowest for the experimentally optimized scales and highest for the BLM method. 
Whereas for the experimentally optimized scales the fit values for the individual a s (M|) 
measurements are perfectly consistent, the consistency for the ECH and the PMS meth- 
ods is only moderate. The results for the ECH and the PMS methods are very similar, 
the correlation p between ECH and PMS scales is almost 1. In the case of the BLM 
method the x 2 / n df indicates that the individual a s (M|) measurements are inconsistent. 
Moreover, the fits using the scales predicted by the BLM method did not converge at all 
for the observables JCEF, O, p D and D$ eneva . Fi gure O shows the correlation between 
the logarithms of the experimentally optimized scales and the logarithms of the scales 
predicted by the ECH, PMS and the BLM methods. For the ECH and the PMS method 
there are significant correlations of p = 0.75 ±0.11. In the case of BLM there is a slightly 
negative correlation of p = —0.34 ± 0.25, compatible with zero within 1.5 a. Our results 
indicate that the ECH and the PMS methods are useful in the case where an experimen- 
tal optimization can not be performed, whereas the BLM method does not seem to be 
suitable for the determination of a s (M|). 
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Observable 


M V P/ 71 IT 1 ) \ 
<*8 ( M z) 


®s (M 2 Z ) 


hit ' H / 71 /T 1 ) \ 

«; t (M 2 Z ) 


M f 7i ir r >\ 


EEC 


0.1142 


0.1133 


0.1135 


0.1142 


AEEC 


0.1150 


0.1063 


0.1064 


0.1179 


JCEF 


0.1169 


0.1168 


0.1169 




1 T 
1 — 1 


U. 110Z 


U. 11U1 


n 1 1 1 1 

U. 1 1 1 1 


U. 1 lOO 





0.1171 


0.1128 


0.1124 






nil 


n 1 1 1 q 

U.lllj 


n 1 194 

yj . x iZ/i 


n 1 144 


BMax 


0.1215 


0.1222 


0.1217 


0.1268 


Bsum 


0.1138 


0.1023 


0.1021 


0.1118 


PH 


0.1215 


0.1197 


0.1198 


0.1258 


Ps 


0.1161 


0.1154 


0.1149 


0.1169 


PB 


0.1172 


0.1190 


0.1203 




u 2 


0.1165 


0.1145 


0.1142 


0.1143 


u 2 


0.1210 


0.1204 


0.1202 


0.1232 


U 2 


0.1187 


0.1110 


0.1108 


0.1118 


T\ Jade 
u 2 


0.1169 


0.1137 


0.1134 


0.1137 


T^Durham 
u 2 


0.1169 


0.1162 


0.1159 


0.1241 


T) Geneva 
u 2 


0.1178 


0.1064 


0.1171 




TpvCambridgc 
U 2 


0.1164 


0.1164 


0.1163 


0.1124 


w. average 
X 2 /n d f 


0.1168 ±0.0026 
6.2 / 17 


0.1154 ±0.0045 
19 / 17 


0.1155 ±0.0044 
19 / 17 


0.1174 ±0.0068 
29 / 13 



Table 4: Comparison of the a s (M|) values obtained using the different methods for 
evaluating the renormalization scale suggested by theory. For each observable the a s (M|) 
values using experimentally optimized scales and a s (M|) values for the scales predicted by 
the PMS, ECH and BLM methods are shown. The errors for the a s (M|) measurements 
are assumed to be identical for all methods (see table ||). The weighted averages are 
calculated using p e s = 0.635 and scaling the errors to yield x 2 / n df = 1 in the case of the 
PMS, ECH and the BLM methods (see text). The x 2 given for the averaging correspond 
to the values before adjusting p e g and rescaling the measurement uncertainties. The fits 
using the scales predicted by BLM did not converge for the observables JCEF, O, pn and 

piGeneva 

u 2 




Figure 11: Comparison between the logarithms of the experimentally optimized renormalization scales and the logarithms of scales 
predicted by PMS, ECH and the BLM method. The vertical error bars indicated represent the uncertainties from the 2-parameter fits 
in X/j, and a s . The vertical bars represent the range of renormalization scale values for the theoretically motivated scale setting methods 
evaluated from the individual bins within the fit range of each distribution, whereas the central value has been derived by considering 
the full theoretical prediction within the fit range. 
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6 Pade Approximation 

An approach for estimating higher-order contributions to a perturbative QCD series 
is based on Pade Approximations. The Pade Approximant [N/M] to the series 

S = S + S lX + S 2 x 2 + ... + S n x n (37) 

is defined || by 

jN/M] = ao + * lX+ u a2X 2 + -\ aNX M ; N + M = n (38) 

L ' J 1 + hx + b 2 x 2 + . . . + b M x M v ; 

and 

[N/M] = S + 0(x N+M+1 ) . (39) 



The set of equations (|39|) can be solved, and by consideration of the terms of 
0(x N+M+1 ) one can obtain an estimate of the next order term Sjv+m+i of the origi- 
nal series. This is called the PA method. Furthermore, for an asymptotic series [N/M] 
can be taken as an estimate of the sum (PS) of the series to all orders. The PA method 
has been used successfully to estimate coefficients in statistical physics |fB3| , and various 



quantum field theories including QCD [pLl|| . Justifications for some of these successes 
have been found in mathematical theorems on the convergence and renormalization scale 
invariance of PAs ||41|| . In many cases the PAs yield predictions for the higher order coef- 



ficients in perturbative series with high accuracy, whereas this accuracy is not expected 
for the lower order predictions such as 0(a\) . For the application of the a s (M|) deter- 
mination from event shapes, the Pade Approximation can serve as a reasonable estimate 
of the errors due to higher order corrections 



For each bin of our observables an estimate for the O(o%) coefficient C(y) can be 
derived from [0/1] with a® = A, b% — —B/A: 



My) 



It should be noted that the PA predictions C Pade (y) are positive by construction which 
will result in large errors for kinematical regions where the 0(a^) contribution is negative. 
The fit range has therefore been determined in the following way: Starting from the same 
fit range as in 0(a 2 s ) , the fit has been accepted if x 2 / n df < 5. Otherwise the fit range 
was reduced bin by bin until the fit yielded x 2 / n df < 5. 

In addition to the 0(a^) fits in the Pade Approximation, the PS method has been used 
as an estimate of the sum of the perturbative series and a s (M|) has been extracted by 
fitting the [0/1] approximation directly to the data. Here, the fit range has been chosen 
to be the same as for the fits in the 0{oP s ) Pade approximation. The \ 2 dependence of 
the a s (M|) fits applying the Pade Approximation as a function of the renormalization 
scale value x M is quite small, especially for the PS method. For most of the observables 
a s (M|) and x^ could not be determined in a simultaneous fit. Therefore, the fits have 
been done choosing a fixed renormalization scale value x M = 1. The uncertainty due to 
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the scale dependence of a s (M|) has been estimated by varying x^ between 0.5 and 2. 

The fit results for the individual observables are listed in tables [5| and ^|. The fit 
applying PS to the Bs U m distribution did not converge for any fit range chosen. For 
the D% eneva distribution, the fits did not converge for either method. Comparing the 
fit results of the O(a^) fits in Pade Approximation with the fit results in 0(a 2 ) applying 

= 1, as given in table |3| , the scale dependence of a s (M|) is reduced for most of the ob- 
servables, as one would expect from measurements using exact calculations in 0(a^) . For 
the PS method, the reduction of the scale dependence is even larger. Here, a s (M|) is less 
scale dependent than in the 0(a 2 ) fits for all observables considered. Figure [12] shows the 
scale dependence of a s (M|) applying the different QCD predictions to the distribution 
of the Jet Cone Energy Fraction as an example. There is almost no \ 2 dependence of the 
a s (M|) fits as a function of the renormalization scale for the PS prediction. For the fits 
applying 0(a^) in the Pade Approximation, the x 2 dependence is less than for the 0(a 2 ) 
prediction. However, the JCEF is one of the few observables, where a simultaneous fit of 
a s (M|) and x M is possible. 



Observable 


Fit Range 


a s 


,(M|) 


Aa s (Scale. 


) Aa s (Tot.) 


EEC 


28.8° - 


151.2° 


0. 


1189 


±0.0016 


±0.0026 


AEEC 


25.2° - 


64.8° 


0. 


.1074 


±0.0030 


±0.0056 


JCEF 


104.4° - 


169.2° 


0. 


.1169 


±0.0006 


±0.0016 


1 -T 


0.07- 


0.30 


0. 


.1207 


±0.0023 


±0.0036 





0.24 - 


0.32 


0. 


.1098 


±0.0014 


±0.0044 


c 


0.32 - 


0.72 


0. 


.1208 


±0.0023 


±0.0039 


BMax 


0.10 - 


0.24 


0. 


.1183 


±0.0016 


±0.0042 


Bsum 


0.14 - 


0.18 


0. 


.1127 


±0.0016 


±0.0068 


Ph 


0.03 - 


0.14 


0. 


.1230 


±0.0015 


±0.0036 


Ps 


0.10 - 


0.30 


0. 


.1252 


±0.0024 


±0.0038 


Pd 


0.07- 


0.30 


0. 


,1045 


±0.0015 


±0.0040 


u 2 


0.05 - 


0.18 


0. 


.1159 


±0.0014 


±0.0042 


d po 

u 2 


0.05 - 


0.18 


0. 


.1199 


±0.0011 


±0.0034 


u 2 


0.10 - 


0.20 


0. 


.1128 


±0.0008 


±0.0030 


D Jade 


0.06 - 


0.25 


0. 


.1142 


±0.0014 


±0.0032 


T~)Durham 
u 2 


0.015 - 


• 0.16 


0. 


.1170 


±0.0009 


±0.0023 


Y) Cambridge 
U 2 


0.011 - 


• 0.18 


0. 


.1164 


±0.0007 


±0.0026 


average 




0.116; 


3± 0.0054 


X 2 /n d f 


= 30/16 



Table 5: Results on a s (M|) for QCD-Fits including the 0(a^) Term in the Pade Ap- 
proximation (PA). For each of the observables the fit range, a s (M|) , the uncertainty due 
to scale variation between 0.5 < x^ < 2 and the total uncertainty are shown. The exper- 
imental errors and the uncertainties due to the hadronization corrections are assumed to 
be the same as for the 0{a 2 s ) measurements. The weighted average is calculated using 
Pes = 0.635 and scaling the errors to yield x 2 / n df = 1 (see text). The \ 2 given for the 
averaging corresponds to the value before adjusting p e fr and rescaling the measurement 
uncertainties. The fit for jjGeneva ^-j no ^. conV erge. 
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Assuming the same correlation as in the 0(a 2 ) fits, the weighted averages of a s (M|) 
over the observables used have been calculated as: 

a B (M|) = 0.1168 ±0.0054 

for the PA fits and 

a B (M|) = 0.1157 ±0.0037 

for the PS fits. The averages are in excellent agreement with the 0(a 2 ) value of 
a s (M|) = 0.1168 ± 0.0026 using optimized scales. The scatter between the observables, 
however, is somewhat larger than in the 0(a 2 ) case. The x 2 / n df for the average values 
is 30/16 and 17/15 respectively. 



Observable 


a s (M'i) 


Aa s (Scale.) 


Aa s (Tot.) 


EEC 


0. 


1147 


±0.0003 


±0.0021 


AEEC 


0. 


1070 


±0.0002 


±0.0048 


JCEF 


0. 


1169 


±0.0003 


±0.0015 


1 - T 


0. 


1165 


±0.0003 


±0.0028 


O 


0. 


1135 


±0.0003 


±0.0042 


c 


0. 


1150 


±0.0003 


±0.0032 


BMax 


0. 


1196 


±0.0003 


±0.0039 


Ph 


0. 


1219 


±0.0004 


±0.0033 


Ps 


0. 


1161 


±0.0003 


±0.0029 


Pd 


0. 


1098 


±0.0003 


±0.0037 


u 2 


0. 


1136 


±0.0003 


±0.0040 


d po 

LJ 2 


0. 


1198 


±0.0003 


±0.0032 


u 2 


0. 


1124 


±0.0003 


±0.0029 


T^Jadc 
u 2 


0. 


1123 


±0.0003 


±0.0029 


T~\Durham 
u 2 


0. 


1164 


±0.0003 


±0.0021 


p. Cambridge 
L> 2 


0. 


1162 


±0.0003 


±0.0025 


average 0.1157 ± 0.0037 


X 2 /n df = 


17 / 15 



Table 6: Results on a s (M|) for QCD-Fits applying the Pade Sum Approximation (PS). 
For each of the observables a s (M|) , the uncertainty due to scale variation between 
0.5 < < 2 and the total uncertainty are shown. The experimental errors and the 
uncertainties due to the hadronization corrections are assumed to be the same as for the 
0(a 2 ) measurements. The weighted average is calculated using p cS = 0.635 and scaling 
the errors to yield x 2 / n df = 1 (see text). The \ 2 given for the averaging corresponds to 
the value before adjusting p e s and rescaling the measurement uncertainties. The fits for 
B S um and ]j^ eneva did not converge. 
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Figure 12: a s (M|) and Ax 2 = x 2 — Xmin f° r the distribution of the Jet Cone Energy 
Fraction as a function of x M from QCD fits applying 0(a 2 ) prediction, 0(al) in Pade 
Approximation and the Pade Sum Approximation. Additionally, the x 2 minimum for the 
0(a 2 ) fit and the renormalization scale value x^ — l have been indicated in the a s (M|) 
curves. 
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7 QCD in the Next-to-Leading Log Approximation 

All orders resummed QCD calculations in the Next-to-Leading Log Approximation 
(NLLA) matched with 0(a 2 ) calculations have been used widely to measure a s (M|) 
from event shape observables 



For a generic event shape observable Y for which the theoretical prediction can be 
exponentiated, the cumulative cross-section at the scale Q 2 = s is defined by: 

Y 

R ^ = \ I % iY ' (41) 



o 



and can be expanded in the form 



R(Y, a s ) = C{a s ) exp E(a s , L) + F(a s , Y) , (42) 
where L = ln(l/Y) and 

oo 

C(a s ) = l + ^CX (43) 

i=l 

OO 2+1 

i=i j=i 



oo 



F(a s ,Y) = Y J F t iY)^ s ■ (45) 

i=l 

The Ci are constant and the Fi(Y) vanish in the infrared limit Y — > 0. The factor E to 
be exponentiated can be written 

S(a s , L) = Lf LL (a s L) + f NL L(a s L) + subleading terms , (46) 

where and /nil represent the leading and the next-to-leading logarithms. They 
have been calculated for a number of observables, including 1 — T f43|, C [Q, B max H45|, 



B sum ^, p H H and D^ urham @, where the NLLA predictions for B max and B s 



entering into this analysis are the recently improved calculations by Yu. L. Dokshitzer et 



al. 4 



Pure NLLA calculations can be used to measure a s (M|) in a limited kinematic region 
close to the infrared limit, where L becomes large. In order to achieve a prediction where 
the kinematical range can be extended towards the 3 jet region, several procedures have 
been suggested to match the NLLA calculations with the calculations in C(a 2 ) . 

The 0(a 2 s ) QCD formula can be written in the integrated form: 

R G(*i) {Y,a s ) = l + A{Y)a s + B(Y )a 2 s , (47) 

where A (Y) and B (Y) are the cumulative forms of A(Y, cos$t) and B(Y, cosi?t) m Eq. 
(EST) , integrated over ■ Together with the first and second order part of Eq. (HI) 



gi{L) = G 12 L 2 + G 11 L 



(48) 



3G 



g 2 (L) = G 23 L 3 + G 22 L 2 + G 21 L , (49) 
the In R matching scheme can be defined as: 

\nR(Y,a s ) = X{a s ,L)+H 1 (Y)a s + H 2 (Y)a 2 s , (50) 

where 

H 1 (Y)=A(Y)-g 1 (L) (51) 

H 2 (Y)=B(Y)- 1 -A 2 (Y)-g 2 (L). (52) 

When combining 0(a 2 ) predictions with NLLA calculations one has to take into ac- 
count that the resummed terms do not vanish at the upper kinematic limit Y max of the 
event shape distributions. In order to correct for this, the resummed logarithms are 
redefined |5tJ by: 

L = \n{l/Y -l/Y max + l). (53) 

Several other matching schemes can be defined which differ in the treatment of the 
subleading terms, thus introducing a principal ambiguity in the matching procedure. The 
In R matching scheme has become the preferred one, because it includes the C 2 and the 
G21 coefficients implicitly and uses only those NLLA terms which are known analytically. 
It yields the best description of the data in terms of x 2 / n df in most cases ||. Therefore, 
we quote the results for the matched predictions using the lni? matching scheme and 
use the R and R — G 2 i matching schemes as defined e.g. in || for the estimation of the 
uncertainty due to the matching ambiguity. 



7.1 Measurement of a s (ilf|) using pure NLLA predictions 

To measure a s (M|) from pure NLLA calculations the fit range has to be restricted to 
the extreme 2-jet region, where L becomes large and the resummed logarithms dominate. 
We define us as the ratio of the resummed logarithms to the non-exponentiating second 
order contributions as follows : 

= £(a 8 ,L) 

H 1 {Y)aL B + H 2 (Y)& a 1 ] 

In addition to the fit range criteria listed in Section ^ we require the minimum of the 
ratio lo over the fit range not to fall below 5 for the fits in pure NLLA. This leads to 
the fit ranges listed in table [7]. For the observable D^ urham the ratio u remains small 
even for small values of y cut . No fit range can be found where the resummed logarithms 
dominate the prediction. Therefore D^ urham has not been used for the fits in pure NLLA. 

Contrary to the 0(a 2 ) predictions an optimization of the renormalization scale (or 
more precisely an optimization of the renormalization scheme) cannot be easily per- 
formed for the resummed NLLA predictions 0]. Therefore the scale was fixed to = 1. 
The uncertainty due to the scale dependence of a s (M|) was estimated by varying the 
scale x M between 0.5 and 2. 
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Observable Fit Range (NLLA) Fit Range (matched) 



1 - T 


0.04 - 


0.09 


0.04- 


0.30 


C 


0.08 - 


0.16 


0.08 - 


0.72 


Rmax 


0.02 - 


0.04 


0.02 - 


0.24 


Bsum 


0.06 - 


0.08 


0.06 - 


0.24 


Ph 


0.03 - 


0.06 


0.03 - 


0.30 


t~\ Durham 
u 2 






0.015 - 


• 0.16 



Table 7: Fit range for the observables in pure NLLA and matched NLLA fits. The 
observable D^ urham has not been used for pure NLLA fits, since no fit range can be 
found, where the resummed logarithms dominate the predictions (see text). 

The fit results for the individual observables are listed in table ||. The weighted aver- 
age of a s (M|) for the 5 observables is 

a s (Ml) = 0.116 ±0.006 

which is in excellent agreement with the average value for the 0(a 2 s ) fits of ot a (M|) = 
0.1168 ±0.0026. 



Observable 




Aq s (exp.) 


Aa s (had.) 


Aa s (seal.) 


Aa s (tot.) 


X 2 / n df 


1 -T 


0.120 


±0.001 


±0.004 


±0.004 


±0.006 


0.59 


C 


0.116 


±0.002 


±0.003 


±0.004 


±0.006 


0.53 


B ma x 


0.111 


±0.004 


±0.003 


±0.002 


±0.006 


2.37 


Bsum 


0.116 


±0.003 


±0.004 


±0.002 


±0.006 


1.24 


Ph 


0.117 


±0.004 


±0.006 


±0.004 


±0.009 


0.43 


average 


0.116 ±0.006 


X z /n 


* = 1-2/4 


Pcff = 


0.71 





Table 8: Results for the a a (M|) fits in pure NLLA for the individual observables together 
with the individual sources of uncertainties and the x 2 / n df for the NLLA fits. The total 
error on a s (M|) listed, is the quadratic sum of the experimental error (statistical and 
systematic uncertainty), the uncertainty due to the hadronization correction and the 
uncertainty due to the scale dependence of a s (M|) . Also listed is the weighted average 
of a s (M|) for the 5 observables together with the x 2 / n df for the averaging procedure 
and the correlation parameter p cS . The x 2 corresponds to the value before readjusting 
according to Eq. [33|. 
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7.2 Measurement of cy. s (M^) using NLLA predictions matched 
with (D(a 2 s ) 

For the QCD fits using NLLA theory matched with 0(a 2 ) predictions the fit range has 
been chosen as the combined fit range for the pure NLLA and the 0{a 2 s ) fits. The results 
for the individual observables in the lni? matching scheme are listed in table |9|. The 
additional uncertainty due to the matching ambiguity has been estimated as the maxi- 
mum difference of a s (M|) in the lni? matching scheme and the two alternative matching 
schemes, R and R — G%\. The dependence on the choice of Y max , i.e. the value of the 
upper kinematic limit for the shape observables used for the redefinition of the resummed 
logarithms according to Eq. Q has been studied by repeating the fits with the value of 
Ymax being reduced by 10 %. The resulting variation in o s (M|) has been found to be 
small, the maximum change is about 1 %. 

The average value of a s (M|) in the lni? matching scheme is 

a s (M|) = 0.119 ±0.005 

which is in good agreement with the average value for the 0(a 2 ) fits of a s (Mf) = 
0.1168 ±0.0026. 

Looking at the individual fit results, one finds from the x 2 / n df that most of the shape 
distributions cannot be successfully described in a fit range expected to apply for the 
combined theory. The a s (M|) values are higher than for the fits in pure NLLA for all 
observables considered. In the case of 1 — T, C and B sum the measured a s (M|) values 
are even above the values for both the pure NLLA fits and the 0(a 2 ) fits using exper- 
imentally optimized scales, where one naively might expect the matched predictions to 
be a kind of 'average' of the individual theories. 

In order to investigate this result further it is instructive to compare the theoretical 
predictions of the shape distributions for the different methods with the data distribu- 



Observable 




Aa s (exp 


.) Aa s (had.) 


Aa s (seal. 


) Aa s (mat.) 


i Aa s (tot.) x 2 /n d f 


1 - T 


0.124 


±0.002 


±0.003 


±0.004 


±0.003 


±0.007 


9.5 


C 


0.120 


±0.002 


±0.002 


±0.004 


±0.004 


±0.007 


15.2 


B ma x 


0.113 


±0.002 


±0.002 


±0.003 


±0.003 


±0.005 


8.4 




0.122 


±0.002 


±0.003 


±0.004 


±0.005 


±0.008 


11.9 


Ph 


0.119 


±0.002 


±0.002 


±0.003 


±0.005 


±0.007 


1.33 


T^Durham 


0.121 


±0.001 


±0.002 


±0.002 


±0.005 


±0.006 


1.70 


average 0.119 ±0.005 


X 2 /n df = 2.3/5 


p cS = 0.57 







Table 9: Results of the QCD fits in the lni? matching scheme for the individual ob- 
servables together with the individual sources of uncertainties and the x 2 / n df for the 
a s (M|) fits. The total error is the quadratic sum of the experimental error (statistical 
and systematic uncertainty), the uncertainty due to the hadronization correction, the un- 
certainty due to the scale dependence of a s (M|) and the uncertainty due to the matching 
ambiguity. Also listed is the weighted average of a s (M|) for the 6 observables together 
with the x 2 / n df for the averaging procedure and the correlation parameter p e g. The x 2 
corresponds to the value before readjusting according to Eq. [33|. 
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tions. Figure [13| shows experimental distributions for 1 — T and C in comparison with 
the fitted curves for three different types of QCD fits, namely 0{a 2 s ) using experimentally 
optimized scales, 0(a 2 ) using a fixed renormalization scale = 1 and the fits in the \nR 
matching scheme. For the fits in 0(a 2 ) using experimentally optimized scales, the data 
are described well over the whole fit range. For the fits in 0(a 2 ) using a fixed renormal- 
ization scale and the fits in the lni? matching scheme, we find only a poor description 
and the slope of both curves show a similar systematic distortion with respect to the 
data. In the case of 0(a 2 ) applying x M = 1 the distortion arises from the wrong choice 
of the renormalization scale. Since the scale value for the matched predictions is also 
chosen to be x M = 1, the similarity of the curves indicates that the subleading and non- 
logarithmic terms originating from the 0(a 2 ) part of the matched theory and introduced 
using the scale value — 1 dominate the InR predictions. It should be noted that the 
matched theory requires a renormalization scale value of O (1). Unlike the 0(a 2 ) case, 
2 parameter fits in a s (M|) and x M do not converge for most of the observables; for such 
low scale values as in 0(a 2 ) the data can not be described at all in the matched theory. 
It seems that the combination of all orders resummed predictions and terms only known 
in 0(a 2 ) results in a systematic shift in a s (M|) due to the impossibility of choosing an 
appropriate renormalization scale value. Although the average values for the 0(a 2 ) fits, 
the fits in pure NLLA and the fits in the lni? matching scheme are in good agreement, 
the matched results should be considered to be less reliable than those of the 0(a 2 ) and 
pure NLLA analyses due to the systematic deviation of the prediction to most of the data 
distributions (see e.g. Figure O and Table 0). 




Figure 13: left part: Comparison of DELPHI data with three different QCD Fits: i ) 0(a 2 ) using an experimentally optimized 
renormalization scale, ii) 0(a 2 ) using a fixed renormalization scale — 1 and Hi) InR matched NLLA (x M = 1) for the thrust 
Distribution. The lower part shows the relative difference (Fit-Data) /Fit. Whereas the 0{a 2 s ) curve describes the data over the whole 
fit range, the slope of the curves for the fixed scale and lni? matching show a similar systematic distortion with respect to the data. g 
right part: The same for the C-Parameter. Here the distortion is even stronger. 
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8 Heavy Quark Mass Effects 



Studies of the influence of quark mass effects on jet cross-sections |HJ have shown 
that these effects can be important in the study of event shape observables and hence 
in the measurement of the strong coupling. For a natural mixture of quark flavours, 
the influence of the quark masses on the measurement of a s (M|) from event shapes is 



expected to be ~ 1% j52| . In order to derive a high precision value of a s (M|) , the O{of) 
measurement with experimentally optimized scales has been refined by considering the 
influence of b quark mass effects in leading order: 



I d 2 a(Y,cos$ T ) 
a tot dYd cos $t 



a s (/i 2 ) • A(F,costf r )[l + Am] + d s 2 (fi 2 ) ■ B(Y, cos $ T ,x^ (55) 



Observable 




Aa s (Mass) 


Aa s (Tot.) 


EEC 





1145 


± 0.0001 


± 0.0028 


AEEC 





1149 


± 0.0001 


± 0.0111 


JCEF 





1180 


± 0.0007 


± 0.0018 


I — T 





1136 


± 0.0001 


± 0.0036 








1184 


± 0.0006 


± 0.0057 


c 





1156 


± 0.0002 


± 0.0036 


B]VIax 





1223 


± 0.0002 


± 0.0041 


Bsum 





1144 


± 0.0003 


± 0.0053 


Ph 





1216 


± 0.0001 


± 0.0060 


PS 





1160 


± 0.0001 


± 0.0029 


Pd 





1196 


± 0.0008 


± 0.0039 


u 2 





1165 


± 0.0002 


± 0.0044 







1207 


± 0.0002 


± 0.0033 


Dl 





1186 


± 0.0001 


± 0.0046 


D Jade 





1182 


± 0.0005 


± 0.0041 


j^Durham 





1172 


± 0.0003 


± 0.0026 


T\ Geneva 
u 2 





1216 


± 0.0013 


± 0.0310 


^-.Cambridge 
U 2 





1176 


± 0.0006 


± 0.0026 


average 0.1174 ± 0.0026 


X 2 / n df = 


6.60 / 17 p eff = 0.615 



Table 10: Results of the refined 0(a 2 ) measurement of a s (M|) including b quark mass 
effects in leading order. The central value of a s (M|) quoted is the average of a s de- 
rived from applying the b pole mass M b and the b running mass m b [Mz) definition. 
Aa s (Mass) = \a s (M b ) — a s {m b (M z )) \ /2 has been taken as an estimate of the uncer- 
tainty due to quark mass effects of higher orders. The total error displayed is the quadratic 
sum of the experimental uncertainty, the hadronization uncertainty, the scale uncertainty 
and the uncertainty due to quark mass effects. The x 2 given, corresponds to the value 
before readjusting according to Eq. [33|. 
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where 



^2 d 2 a q (Y, cos T )/ ^2 a i 



Am = ^ - 1 (56) 

d 2 a q (Y,cos9 T )/ ^ a q 

g,mi,=0 q,mi,=0 



The coefficients Am for the 18 observables studied have been computed numerically 
for two different definitions of the b-quark mass: The b pole mass of Mb ~ 4.6 GeV/c 2 
and the b running mass at the Mz, m^iMz) ~ 2.8 GeV/c 2 j5T|. All definitions of the 
quark mass are equivalent to leading order; differences are entirely due to higher orders in 
a s . a s (M|) has been determined applying both mass definitions. The average of a s (M|) 
derived from the pole mass definition and the running mass definition has been taken as 
the central value, and half the difference has been taken as an estimate of the uncertainty 
due to higher order mass effects. The results for the individual observables are listed in 
table HO. 



9 Summary 

From 1.4 Million hadronic Z decays recorded with the DELPHI detector and repro- 
cessed with improved analysis software, the distributions of 18 infrared and collinear safe 
observables have been precisely measured at various values of the polar angle $t of the 
thrust axis with respect to the beam direction. The dependence of all detector proper- 
ties has been taken fully into account to achieve the best possible experimental precision. 
In order to compare with QCD calculations in 0(a 2 ) , hadronization corrections are eval- 
uated from precisely tuned fragmentation models. 

The precise data are used to measure a s (M|) applying a number of different methods 
described in the literature. The most detailed studies have been performed in second 
order pertubative QCD. Fits taking explicit account of the a s dependent event orienta- 
tion as predicted by QCD with experimental acceptance corrections less than ~ 25% and 
hadronization corrections less than ~ 40% yield the result that the data can be surpris- 
ingly well described in 0(a 2 ) by using a common value of a s (M|) with a small uncertainty. 

Taking account of the correlation among the observables an average value of a s (M§) = 
0.1168 ± 0.0026 is obtained from the data. The consistency of the individual a s (M|) de- 
terminations from the 18 shape distributions is achieved by using the different values of 
the renormalization scales as obtained from the individual fits, i.e. applying the so called 
experimental optimization method. The significance of the fits is improved due to the 
large number of data points per distribution. Thus, definite results concerning the choice 
of the optimal renormalization scale value become possible. It should be pointed out that 
for most of the investigated observables the scale dependence of a s is very small in the 
vicinity of the experimentally optimized scale. The quoted error of a s (M|) includes the 
uncertainty due to a variation of the experimentally optimized scale in the range between 
0.5 • x^ p and 2. • xf, xp . 

An analysis with a fixed renormalization scale value of — 1 yields an unaccept- 
able description of the data for many observables and leads to a wide spread of the 
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<x,(M|) values. In contrast to fits applying experimentally optimized scales, the stability 
of cc s (M|) with respect to the choice of the fit range is in general quite poor. Due to the 
improved accuracy of the data, systematic differences between the 0(a 2 ) analysis using 
experimentally optimized and fixed renormalization scale values x^ — 1 become visible. 
These differences propagate also into the matched NLLA-analysis. 

To check the reliability of the a s results obtained from the experimentally optimized 
scales three further approaches for choosing an optimized value of the renormalization 
scale have been investigated: The principle of minimal sensitivity (PMS), the method 
of effective charges (ECH), and the method of Brodsky, Lepage and MacKenzie (BLM). 
The weighted averages of a s from the three methods are in excellent agreement with the 
weighted average of a s obtained from the experimental optimization, but their scatter is 
larger. The scatter is largest for BLM. A significant correlation between the renormaliza- 
tion scale values evaluated with ECH and PMS with the experimentally optimized scale 
values is observed. No such correlation exists for the BLM scales. 

Further approaches to estimate the influence of higher order contributions to the per- 
turbative QCD series are based on Pade approximants. The [0/1] Pade approximant has 
been used as an estimate of the sum of the perturbative series as well as for the extrap- 
olation of the unknown O(a^) coefficients for the 18 distributions. In both studies the 
renormalization scale has been set to = 1. Again the average values of a s (M|) are 
consistent with the average value from the experimental scale optimization in 0(a 2 ) . 

While all above mentioned determinations of a s (M|) use fixed order perturbation 
theory the last part of the paper describes measurements of a s (M|) using all orders re- 
summed calculations in the next-to- leading logarithmic approximation (NLLA). In a first 
step, pure NLLA predictions have been confronted with the data in a limited fit range 
where the ratio of the resummed next-to-leading logarithms to the non-exponentiating 
0(a 2 ) contributions is large. The very good agreement between the average value of 
a a (M|) obtained from the pure NLLA fits with a renormalization scale value x^ — 1 
and the 0(a 2 ) fits using experimentally optimized scales is remarkable. In a further step 
NLLA matched to 0(a 2 ) calculations have been applied. The corresponding average 
value of a s (M|) is again consistent with the 0(a 2 ) result though the a s values from all 
investigated observables are systematically higher. More importantly, the application of 
matched NLLA to the high precision data reveals that the trend of the data deviates in 
a systematic fashion from the predictions of the matched theory. This problem has not 
previously been observed. The matched results should be considered to be less reliable 
than those of the 0(a 2 ) and pure NLLA analyses. 

The a s values derived from the different approaches considered are in very good agree- 
ment. The 0(a 2 ) analysis applying a simultaneous fit of a s and x^ to the experimental 
data yields superior results in all respects. 

In a final step the influence of heavy quark mass effects on the measurement of a s (M|) 
has been studied. The weighted average of a s (M|) from the 0(a 2 ) measurements using 
experimentally optimized renormalization scale values and corrected for the b — mass to 
leading order yields 



a s (M 2 z ) = 0.1174 ±0.0026. 
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prediction M^z) x v X 2 / n df 

0(af) exp. opt. scale 0.1169 ± 0.0017 0.0820 1.05 

ECH/FAC 0.1169 ±0.0017 0.2189 2.75 

PMS 0.1168 ±0.0017 0.1576 1.91 

Pade 0(a 3 s ) (fixed scale) 0.1169 ± 0.0016 1.0 3.12 

Pade Sum (fixed scale) 0.1169 ± 0.0015 1.0 3.05 

Pade 0(a 3 s ) (exp. opt. scale) 0.1164 ± 0.0015 0.1814 2.45 

Q{a 2 s ) (fixed scale) 0.1191 ± 0.0024 1.0 7.7 

0(aj) exp. opt. scale 0.1180 ± 0.0018 
+ LO quark mass effects 



Table 11: Summary of a a (M|) measurements from the distribution of the Jet Cone 
Energy Fraction (JCEF). 



The result is consistent with the result of a previous DELPHI publication [fj] if com- 
pared with partonshower hadronization corrections. 

Among the observables studied, the Jet Cone Energy Fraction (JCEF) Jl9| natu- 
rally reveals some superior properties. First, the size of the hadronization correction is 
extremely small, the average correction within the applied fit range being only about 
3.5%. Furthermore, the second order contribution to the cross-section is quite small. 
Within the applied fit range, the average ratio of the second to first order contributions 
is {tnlo) ~ 25% if x M is fixed at 1 and only {r^io) ~ 6% if the experimentally optimized 
scale value is applied. This indicates a good convergence behavior of the corresponding 
perturbative series. Also the scale dependence of a s (M|) derived from JCEF is very 
small. Both experimental and theoretical uncertainties are smallest for the a s (M|) mea- 
surement from the JCEF distribution. 



Table |TT] shows a summary of a s measurements from the JCEF distribution for the 
different methods. The 0(a 2 s ) fit of a s (M|) applying a fixed renormalization scale value 



x 



1 clearly fails to describe the data with x / n df = 965/125. However, even for 
this method, the deviation of a s (M|) from the value obtained using the experimentally 
optimized scale value is only about 2 %. All other approaches yield nearly identical a s 
values within a few per mille. The deviation of the a s values derived from the different 
methods, which can serve as an estimate of the theoretical uncertainty due to missing 



higher order terms |]35| , is clearly smaller than the uncertainty of ±0.0008 derived from 
the variation of the renormalization scale value. Due to the outstanding qualities of this 
observable, the JCEF is considered as best suited for a precise determination of a s (M|). 
After correcting the measured value for heavy quark mass effects, the final result is 



a s (M|) = 0.1180 ± 0.0006(exp.) ± 0.0013(hadr.) ± 0.0008(scale) ± 0.0007(mass). 

Comparing this result with other recent precision measurements of a s (M|) , there is 
very good agreement with the determination of a s (M|) = 0.1174 ± 0.0024 from Lattice 
Gauge Theory []53| and the recent result from an NNLO analysis of ep deep inelastic 
scattering data of a s (M|) = 0.1172 ± 0.0024 p§. The result is also in good agreement 



with the result from the LEP electroweak working group of a s (M|) = 0.119 ± 0.004 [56 
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from the standard model fit to the full set of electroweak precision data. Compared with 
the most recent result from spectral functions in hadronic tau decays, a s is smaller than 
the central value of a s (M|) = 0.1219 ±0.0020 quoted in |]57|, but in very good agreement 
with the value of a s (M|) = 0.1169±0.0017 derived in |57| by using an alternative analysis 
method considering renormalon chains. 
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